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Abstract. We construct some new Integrable Systems (IS) both classical and 
quantum associated with elliptic algebras. Our constructions are partly based 
on the algebraic integrability mechanism given by the existence of commuting 
families in skew fields and partly - on the internal properties of the elliptic 
algebras and their representations. We give some examples to make an evidence 
how these IS are related to previously studied. 



Introduction. This paper is an attempt to establish a direct connection between 
two close subjects of modern Mathematical Physics - the theory of Integrable 
Systems (IS) and the Elliptic Algebras. The aim of this connection is two-fold: 
we clarify some our recent results in the both domains and fill the natural gap 
proving that some large class of the Elliptic Algebras carries in fact the structure 
of an IS. 

We will start with a short account in the subject of the story and will describe 
briefly a type of the IS's under cosideration. 

In [12], B.Enriquez and second author had proposed a construction of com- 
muting families of elements in skew fields. They explained how to use this con- 
struction in Poisson fraction fields to give an another proof of the integrability of 
Beauville - Mukai integrable systems associated with a K3 surface S ([!]). 

The Beauville-Mukai systems had appeared as the Lagrangian fibrations which 
have the form S^^'^ — > |£| = P(iJ°(S', £)), where S"'^' is the Hilbert scheme of g 
points of S, equipped with a symplectic structure introduced in [19], and C is 
a line bundle on S. Later, the authors of [8] explained that these systems are 
natural deformations of the "separated" (in the sense of [13]) versions of Hitchin's 
integrable systems, more precisely, of their description in terms of spectral curves 
(already present in [15]). Beauville-Mukai systems can be generalized to surfaces 
with a Poisson structure (see [3]). When 5* is the canonical cone Cone(C) of an 
algebraic curve C these systems coincide with the separated version of Hitchin's 
systems. The paper [12] shows how the commuting families construction provides 
a quantization of these systems on the canonical cone. 

A quantization of Hitchin's system was proposed in [2]. It seems interesting to 
construct quantization of Beauville-Mukai systems and to compare it with Beilin- 
son - Drinfeld quantization. We had conjectured the correspondence between the 
results of [2] and a quantization of fraction fields in [9] . A part of this program 
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was realized in the [9] for the case of 5" = r*P^, where means a rational curve 
with n marked points. 

Another main object of the paper is the famihes of Elliptic Algebras. These 
algebras (with 4 generators) were appeared in the works of Sklyanin [20], [21] 
and then they were generalized (for any number of generators) and intensively 
studied by Feigin and one of the authors([22],[27],[28]). These algebras can be 
considered as a flat deformation of the polynomial rings such that the linear ( in 
the deformation parameter ) term is given by the quadratic Poisson brackets. The 
geometric meaning of the Poisson structures was underscored in [29] (see also [41] 
and [36]) - it is the Poisson structures of moduli spaces of holomorphic bundles 
on the elliptic curve. We will use the recent survey [30] as a main source of the 
results and references in the theory of eUiptic algebras. 

The relevance of Elliptic Algebras to the theory of IS was manifested from the 
very beginning. They appeared in Sklyanin's demonstration of the integrability 
of the Landau-Lifshitz model ([20],[21])in the frame of the Faddeev's school ide- 
ology (Quantum Inverse Scattering Method and R— matrix approach). Later 
Cherednik had observed the relation of the Elliptic Algebras defined in [22] with 
the Belavin i?— matrix (see [32]). An interesting observation of Krichever and 
Zabrodin giving an interpretation of a generator in the Sklyanin's Elliptic Alge- 
bra as a Hamiltonian of 2-point Ruijsenaars IS ([16]) was generalized later in [4] 
to the case of double-elliptic 2-point classical model. However, all apphcations of 
these algebras to the IS theory had an indirect character until the last two years. 

We should mention also the results which were obtained in the paper of Sokolov- 
Tsyganov ([34]) where they construct, (using the Sklyanin's definition of quadratic 
Poisson structures), some classical commuting families associated with this Pois- 
son algebras. The integrabihty of these families is imphed by the Sklyanin's 
ideology of Separated Variables (SoV) and technically is based on some gener- 
alization of the classical methods going back to Jacobi, Liouville and Stackcl ( 
which ideologically is very close to the classical part of the theorems in [12] and 
[35] ). However, all their results are in the " non-eUiptic" case. 

Recently, an example of integrable systems associated with linear and qua- 
dratic Poisson brackets given by the elliptic Belavin-Drinfeld classical r-matrix 
was proposed in [18]. This system (an elliptic rotator) appears both in finite 
and infinite-dimensional cases. They give an elliptic version of 2D ideal hydro- 
dynamics on the symplectomorphism group of the 2-dimensional torus as well 
as on a non-commutative torus. We should also mention another appearance of 
the elliptic algebras in the context of Non-Commutative geometry (see [40]). It 
would be interesting to relate them to a numerous modern attempts to define a 
Non-Commutative version of IS theory. 

In this paper we construct some IS which appear directly in the frame of the 
elliptic algebras. The elliptic algebras are figured here in two-fold way: sometimes, 
they are carrying the commuting families of Hamiltonians, sometimes, - they 
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provide a necessary background to our constructions which use their properties 
and representation theory ( basically, the so-called "functional realization" and 
the " bosonization" mappings). 

A family of coordinated (compatible) elliptic Poisson structures was introduced 
in [31]. This family contains three quadratic Poisson brackets such that their 
generic linear combination is the quasi-classical limit g„(^) of an elliptic algebra 
Qn{£-i v)- The famous Lenard - Magri scheme provides the existence of a classical 
integrable system associated with the elliptic curve but it was not clear how to 
get a quantum counterpart of this system because of lack of the knowledge how 
to quantize in general the Magri-Lcnard scheme. Nevertheless, this system can 
be quantized for n = 2m using the approach developing in [12] and it is one of 
the main results of our paper. 

Some of the "elliptic" commuting elements which we construct in this arti- 
cle are related to the quantum version of the above-mentioned bi-hamiltonian 
system. Some other families are associated with a special choice of the elliptic 
algebra. These famihes are obtained, in one hand, as the direct apphcation of 
the construction from [12] to the elliptic algebras and, in other hand, by using 
the properties of the "bosonization" homomorphism, constructed in earlier works 
([2 7], [28]). Some of these families (under the appropriate choice of their numeric 
parameters) may be interpreted as examples of algebraic completely integrable 
systems. We give a geometric interpretation to some of them describing a link 
with the Lagrangian fibrations on symmetric products of elliptic curve cone, giv- 
ing a version of the Bcauvillc-Mukai systems (sec [1], [13], [12], [39]). 

Roughly speaking, the integrable systems of the first type have as the phase 
space a 2m— dimensional component of the moduli space of parabolic rank two 
bundles on the given elliptic curve £. More precisely, the coordinate ring of 
the open dense part of the latter has a structure of a quadratic Poisson algebra 
isomorphic to q2m{^)- We have explicitly verified that the quantum commuting 
elements from our construction are the same as the latter obtained from the 
Lenard- Magri scheme for m = 3 (the first non-trivial case). 

Our main theorem (Thm.3.1) takes place for n = 2m, but we are sure that 
there are some interesting integrable quantum systems in the case of n = 2m + 1. 
It would be interesting to study the bi-hamiltonian structures giving the algebra 
Q2m+i{^) using the results of Gelfand-Zakharevich ([33]) on the geometry of bi- 
hamiltonian systems in the case of impair-dimensional Poisson manifolds. The 
precise quantum version of these systems in the context of the elliptic algebras 
Q2m+i{^,v) is still obscure and should be a subject of further studies. 

The theorems from [12] may be also interpreted as an algebraic version of the 
SoV method ( as it was argued in [35]). Hence, it is very plausible that some of 
our quantum commuting families arising from the generalization of the Jacobi- 
Liouville-Stackel conditions ( which are guaranteed by existence conditions of the 
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Cartier - Foata NC determinants) are the quantum elliptic versions of the IS from 
[34]. We hope to return to this question in our future paper. 
We give also some low-dimensional examples of our construction. 

1. Commuting families in some non- commutative algebras 

Let A be an associative algebra with unit. We will suppose further that we 
work with a skew field. 

Fix a natural number n > 2. We will suppose that there are n subalgebras 
Bn C Am A such that for any couple of the indices j,l <i,j <n the elements 
G Bi and G Bj arc commuting (and inside of the each subalgebra Bi the 
elements generically arc not commuting). 

Let us consider the following data: take an n x (n + 1) matrix M. 
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such that all elements of i— th row belong to the i—ih subalgebra 

We will denote hj Ai^ a n x n matrix which is resulted from the matrix Ai 

with the i — th row erased. The corresponding determinants (if they exist) will 

be denote by M\ 

Now we observe that all principal n + 1 minors 

M°,M\... ,M" 

of the highest order n are correctly defined with the help of Cartier - Foata 
non- commutative determinant. 

It's definition repeats verbatim the standard one because in all n matrices 

M^,M\... .M" 

the entries in the different rows arc commute and each summand in the standard 
definition of the determinant contains the product of n elements of different rows 
whose product is order - independent. 
The following theorem was proved in [12] 

Theorem 1.1. Suppose that one of the matrix A4\i = 1,... ,n, say, Ad'^, is 
invertible. Then the elements Hi — {M'^)~^M'^,i — 1, . . . ,n are commute in the 
skew field A (which we are identifying with its image in the fraction field Frac{A) 
under the monomorphic embedding A Frac( A) . 

The proof of the theorem is achieved by some tedious but a straightforward 
induction procedure. 

The similar results were obtained in the framework of so-called multi-parametric 
spectral problems in Operator Analysis ([38]) and in the framework of so called 
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Seiberg-Witten integrable systems associated with a hyperelliptic spectral curves 
in [35]. 

The important step in the demonstration of the 1.1 is the following "triangle" 
relations which are similar to the usual Yang -Baxter relation: 

M^(M°)-^M^' = M^{M^)-^M' (2) 

B'^{M^)-^B^^ = B'\M^y^B'\ (3) 

where B"^^ is the co-factor of the matrix element 6*-, < i,j < n. 
The theorem 1.1 can be re-formulated to give the following result 

Corollary 1.1. Let A be an algebra, {fi,j)o<i<n,i<j<n be elements of A such that 

fi,jfk,£ fkyifij 

for any i, j, k, £ such that j ^ £. For any / C {0, . . . , n}, J C {1, . . . , n} of the 
same cardinality, we set A/^j = So-eBij(/ j) Here Bij(7, J) denotes a 

set of bijections between two sets of the indices I and J. Assume that the A/ j 
are all invertible. Set Aj = '^{i^,,,^n},{o,...,i,...,n} be the element corresponded to 
J with the index i omitted. Then the 

H, = (Ao)-^A, 

all commute together. 

1.1. Poisson commuting families. We will fix a base field k of characteristic 
7^ 2. The following observation is straightforward: 

Lemma 1.1. If B is an integral Poisson algebra, then there is a unique Poisson 
structure on Prac(S) extending the Poisson structure of B. 

This structure is uniquely defined by the relations 

{1//, g] - -{/, 9}/f\ {1//, i/g] = {/, ^}/(/V) 

Theorem 1.1 has a Poisson counterpart. 

Theorem 1.2. Let A be a Poisson algebra. Assume that A is integral, and 
let TT : A ^ Frac(y4) be its injection in its fraction field. For each n—uple of 
Poisson subalgebras Bi, . . . ,Bn of A such that the elements of pair-wise different 
subalgebras Bi are Poisson commuting (for any pair of indices i.j.i ^ j the 
elements bi G Bi and bj G Bj we have {bi,bj} = 0.) We will write the analogue 
of the matrix (1) like the vector-row: = [6*^, 6^, . . . , 6"] , where 

I ^1 \ 

b'^ : . 
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We set 

Af^' = det[6°,... ,6"]. 
Here as usual we denote by U the i—th omitted column. Then if A'^'^^^ is nonzero 
we setHf^^ — Af'^'^Y^o '''''' o^^^ the family {Hf'^'''')i=i^,„^n is Poisson-commutative: 

for any pair (i,i). 

Remark 1. The elements 6f and 6^- of the matrix M. belongs to the different 
subalgebras Bi and Bj ii i ^ j and hence are Poisson-commute. This condition 
reminds the classical constrains on the Poisson brackets between matrix elements 
appeared in XIX century in the papers of Stockel on the Separation of Variables of 
Hamilton- Jacobi systems ([37]). So our theorem can be considered as an algebraic 
re-definition of the Stackel conditions 

1.1.1. Plucker relations. We want to remind the important step of the second 
proof in [12] which shows the relations between the commuting elements and the 
Plucker-like equations. 
We have to prove 

^ciass|^ciass^ ^ciass| ^ ^^^^-^ permutation in {i, j, k) = 0. (4) 

We have 

n n 
p=l a=0 

where (if i) 

^[i...p...n) ^ ^_^^i„<, det[6° . . . 6" . . . 6\ . . 6"](f) 

(which means that the p— th row in the matrix \fp, . . . ,b", . . . , 6' ... 6"] should 
be erased.) 

We set la<i = 1 if q; < i and otherwise. If o; = i we assume Aa,i — 0. Now 
we have 

n n 
p=l a,/3=0 

SO identity (4) is a consequence of 

V(i,j, A;,a,/5,7), ^ e(cT)A«,^(j)A;3,^(j)A^,^(fc) = 0. (5) 

(T6Perm(i,j,fe) 

When card {a, ... , A;} = 3, this identity follows from the antisymmetry relation 
+ = 0. 

When cardja, . . . , fc} = 4 (resp., 5,6), it follows from the following Plucker 
identities (to get (5), one should set V = (yl®")®" and A some partial determi- 
nant). 
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Let V he a vector space. Then 

- if A e A^iV), and a, b,c,d e V, then 



A(a, 6)A(c, d) - A(a, c)A{b, d) + A(a, d)K{h, c) = 0; 



(6) 



- if A e A^{V) and a, 6, c, 6', d e then 

A(6, c, c')A(a, c, 6')A(6, 6', c ) + A(6, c, 6')A(c, h\ c')A(a, 6, c ) 
- A(6, c, 6')A(a, c, c')A(6, h\ c') - A(6, c, c')A(c, 6', c')A(a, 6, 6') = 0; 

- if A e A^(l^) and a, 6, c, a', 6', d e 1/, then 

A(6, c, 6', c')A(a, c, a', c')A(a, 6, a', 6') + A(6, c, a', c')A(a, c, a', 6')A(a, 6, 6', c') 
+ A(6, c, a', 6')A(a, c, 6', c')A(a, 6, a', d) 

- A(6, c, 6', c')A(a, c, a', 6')A(a, 6, a', c') - A(6, c, a', 6')A(a, c, a', c')A(a, 6, 6', c') 

- A(6, c, a', c')A(a, c, 6', c')A(a, 6, a', 6') = 0. (7) 

We refer to [12] for the proof of these identities. We will need them below in 
some special situation arising with the commuting elements in associative and 
Poisson algebras which are directly connected with elliptic curves and vector 
bundles on them. This Plucker relations can be interpreted as a kind of Riemann- 
Fay identities which are in its turn related to integrable (difference) equations in 
Hirota bilinear form. 



Now we should describe one of the main heroes of our story - the elliptic 
algebras. We will follow to the survey [30] in our notations and also we will refer 
to it as a main source of the results and its proofs in this section. 

2.1. Definition and the main properties. The elliptic algebras are the asso- 
ciative quadratic algebras Qn,ki^, v) which were introduced in the papers [22, 28]. 
Here S is an elliptic curve and n, k are integer numbers without common divisors 
,such that 1 < k < n while r/ is a complex number and Qn,k{^, 0) = C[a;i, ...,Xn]- 
Let S = C/r be an elliptic curve defined by a lattice F = Z©tZ, t e C, '^t > 0. 
The algebra Qn,k{^,v) has generators Xi,i e Z/nZ subjected to the relations 



and have the following properties: 1) Qn,k{^,'>]) = C © Qi © Q2 © ••• such that 
Qa * Q/3 = Qa+i3, here * denotes the algebra multiplication. In other words, the 
algebras Qn^k{^,r]) are Z - graded; 
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2) The Hilbert function of Qn,k{^,v) is ^^>q dim Qa'^" = (i-t)n ■ 

Wc consider here the theta-functions 9i{z),i = 1, . . . , n as a base in the space of 

the theta-functions 6„(r) of the order n which are subordinated to the following 

relations of quasi-periodicity 

ei{z + l) = e,{z), e,{z + T) = {-l)^exp{-2TTy^{-l)nz)9i{z),i = 0,... ,n-l. 

The theta-function of the order 1 6{z) G 0i(r) satisfies to the conditions 
^(0) = and 9{-z) =9{z + t) = - exp(-27r/( - l)z)9{z). 

We see that the algebra Qn,k{^,rj) for fixed £^ is a flat deformation of the 
polynomial ring C[a;i, x^]- The linear (in ?7)term of this deformation gives rise 
to a quadratic Poisson algebra qn,k{^)- 

The geometric meaning of the algebras Qn,k was underscored in [29], [41] where 
it was shown that the quadratic Poisson structure in the algebras qn,k{^) associ- 
ated with the above-mentioned deformation is nothing but the Poisson structure 
on P""^ = ¥Ext^{E, O), where E is a, stable vector bundle of rank k and degree 
n on the elliptic curve £. 

In what follows we will denote the algebras Qn,i{£, r]) by Qn{£: v) 

2.2. Algebra Qn{S,r])- 

2.2.1. Construction. For any n G N, any elliptic curve S = C/P, and any point 
77 G ^ we construct a graded associative algebra Qn{S, v) = C © Fi © F2 © . . . , 
where Fi = e„(P) and = -S"e„(P). By construction, dimF^ = n(n+i)..jn+a-i) _ 
It is clear that the space can be realized as the space of holomorphic symmetric 
functions of a variables {f{zi,..., Za)} such that 

f{zi + 1,Z2, . . . , Za) — f{zi, . . . ,Za), 

f{z, + T,Z2,..., 2„) = (-l)"e-^™/(^i, ...,Za). 



(8) 



For f E Fa and g E Fp we define the symmetric function f*g of a + p variables 
by the formula 

f*g{zi, Za+p) = ^ ^ f{z^^+Pr], ... , z^^+Pr))g{z^^^^-ar], ... , z^^^^-ar))x 

X n ''XXT - 

a + l < j < a + P 

In particular, for f,gEFi we have 
. , X , X , .9{zi - Z2 - nr]) 9{z2- Zi-n-q) 

9[Zi-Z2) 9[Z2-Zi) 

Here 9{z) is a theta function of order one. 
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Proposition 2.1. If f & and g G F^, then f * g E Fa+/3. The operation * 

defines an associative multiplication on the space (Ba>oFa 

2.2.2. Main properties of the algebra Qn{S,ri). By construction, the dimensions 
of tiie graded components of tlie algebra Qn{S,r]) coincide witli tliose for tlie 
polynomial ring in n variables. For 77 = the formula for f * g becomes 

f * g{Zi, . . . ,Za+l) ^ -^^^ ^ f{Zai, . . . ,Za„)g{Za„_^_^, . . . ,Za^_^_^). 

This is the formula for the ordinary product in the algebra S'*0„(r), that is, 
in the polynomial ring in n variables. Therefore, for a fixed elliptic curve £ 
(that is, for a fixed modular parameter r) the family of algebras Qn{S,ri) is a 
deformation of the polynomial ring. In particular, there is a Poisson algebra, 
which we denote by One can readily obtain the formula for the Poisson 

bracket on the polynomial ring from the formula for f * g by expanding the 
difference f*g — g*f in the Taylor series with respect to r;. It follows from the 
semi-continuity arguments that the algebra Qn{£, v) with generic 77 is determined 
by n generators and quadratic relations. One can prove (see §2.6 in [30]) 

that this is the case if 77 is not a point of finite order on £, that is, Nrj ^ F for 
any iV G N. 

The space On(r) of the generators of the algebra Qn{S,ri) is endowed with 
an action of a finite group r„ which is a central extension of the group T/nF of 
points of order n on the curve £. It immediately follows from the formula for the 
product * that the corresponding transformations of the space Fa — <S'"©„(r) are 
automorphisms of the algebra Qn{S,r]). 

2.2.3. Bosonization of the algebra Qn{S,r]). The main approach to obtain rep- 
resentations of the algebra Qn{S,r]) is to construct homomorphisms from this 
algebra to other algebras with simple structure (close to Weil algebras) which 
have a natural set of representations. These homomorphisms are referred to as 
bosonizations, by analogy with the known constructions of quantum field theory. 

Let Bp^niv) be a Z^'-gradcd algebra whose space of degree (cti, . . . , Up) is of the 
form {f{ui, . . . ,'Up)e"^ . . ■Cp''}- where / ranges over the mcromorphic functions 
of p variables and elements of the algebra Bp^n{v)- Let Bpn{r]) 

be generated by the space of meromorphic functions f{ui, . . . ,Up) and by the 
elements ei, . . . , ep with the defining relations 

Bafiui, ...,Up) = f{ui -2r],...,Ua + {n- 2)ri, ...,Up- 2ri)ea, 
ea^is = e/3e«, f{ui, Up)g{ui, ...,Up)= g{ui, Up)f{ui, ...,Up) ^^q-^ 

We note that the subalgebra of Bp^niv) consisting of the elements of degree 
(0, . . . , 0) is the commutative algebra of all meromorphic functions of p variables 
with the ordinary multiplication. 
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Proposition 2.2. Let rj G S be a point of infinite order. For any p G N there is a 
homomorphism 0^: Qn{S,ri) — >• Bp^niv) that acts on the generators of the algebra 
Qn{S,ri) by the formula : 



Here f G 0„(r) is a generator of Qn{S,r]) and the product in the denominator is 
of the form Yl-_^^ 9{ua - Ui). 

2.2.4. Symplectic leaves. We recall that Qn{^, 0) is the polynomial ring 5'*0„(r). 
For a fixed eUiptic curve E — C/F we obtain the family of algebras Qn{£,ri), 
which is a fiat deformation of the polynomial ring. We denote the corresponding 
Poisson algebra by qn{S). We obtain a family of Poisson algebras, depending on 
£, that is, on the modular parameter r. Let us study the symplectic leaves of 
this algebra. To this end, we note that, when passing to the limit as r] ^ 0, 
the homomorphism (pp of associative algebras gives a homomorphism of Poisson 
algebras. Namely, let us denote by bp^n the Poisson algebra formed by the ele- 
ments Eai,...,aj,>o fau-,apiui, ... , Up)e"' . . . Cp", whcrc /ai,...,ap are meromorphic 
functions and the Poisson bracket is 

{ua,Uf3} = {ea,ep} = 0; {ea-,up} = -2ea; {ca, = (n - 2)eo:, 

where a ^ [3. 

The following assertion results from Proposition 6 in the limit as r/ — > 0. 

Proposition 2.3. There is a Poisson algebra homomorphism i/jp'. qn{S) bp^n 
given by the following formula: if f & 0„(r), then 



Let {9i{u);i G Z/nZ} be a basis of the space 0n(r) and let {xi;i G Z/nZ} 
be the corresponding basis in the space of elements of degree one in the algebra 
Qn{£,r]) (this space is isomorphic to 0„(F)). For an eUiptic curve £ C P*~^ 
embedded by means of theta functions of order n (this is the set of points with 
the coordinates (^^o(^) : • • • : 9n-i{z))) wc denote by Cp£ the variety of p-chords, 
that is, the union of projective spaces of dimension 1 passing through p points 
of £. Let K{Cp£) be the corresponding homogeneous manifold in C". It is clear 
that K{Cp£) consists of the points with the coordinates 




(11) 





l<a<p 



, where Bq, G C. 
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Let 2p < n. Then one can show that dim K{CpS) = 2p and K{Cp^i£) is the 
manifold of singularities of K{Cp£). It follows from Proposition 7 and from the 
fact that the Poisson bracket is non-degenerate on 6p „ for 2p < n and Cq, 7^ that 
the non-singular part of the manifold K{CpE) is a 2p- dimensional symplectic leaf 
of the Poisson algebra qn{S)- 

Let n be odd. One can show that the equation defining the manifold K{Cn-i£) 
is of the form C = 0, where C is a homogeneous polynomial of degree n in the 
variables Xj. This polynomial is a central function of the algebra qn{S)- 

Let n be even. The manifold K{Cn^£) is defined by equations Ci — and 
C2 = 0, where degCi = degC2 — n/2. The polynomials Ci and C2 are central in 
the algebra qn{£)- 



There are (at least two) different ways how to construct some commuting ele- 
ments and IS associated with the elliptic algebras. We will start with the general 
statements about the commuting elements arising from the ideas and construc- 
tions of the section 1 . 

3.1. Commuting elements in the algebras Qn{£]r}). Let us consider the 
following Wcyl-likc algebra V„ with the set of generators fi, ■ ■ ■ , fn, Zi, . . . ,Zn 
subjected to the relations 

= [fi, fj] = [zi, Zj] = [fi, Zj] {i ^ j), fiZi = {zi - nri)fi. 

We assume the following commutation relations between the functions in the 
variables Zi and the elements fj: 

fjF{zi, ... , z„) = F{zi, . . . ,Zj -nrj,... , Zn)fj. 

We should precise that the algebra V„ is spanned as a vector space by the 
elements of the form F[zi, . . . , Zn)f^'^ . . . /™", where F is a meromorphic function 
in n variables. 

Remark 2. We should observe also that the algebra V„ looks different from the 
above-mentioned Weyl-like algebras Bp^n but it is isomorphic to the algebra Bn,,n 
and may be reduced to it by a change of the generators. We will return below to 
a geometric interpretation of the algebra V„. 

Now we take the following determinant 



3. INTEGRABLE SYSTEMS 



^o(^i) ^^l(^l) 

^0(^22) ^1(^22) 



On-2{.Zi) 
On-2{z2) 



^n-l(^l) 
On-l{z2) 



Oo{Zn) Oi{Zn) 



^n-2{Zn) 
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cexjp{z2 + 2zs + ... + {n- 1)^„) Y\_ ^i^i ~ ^i)^ 

l<i<j<n i=l 

where the constant c is irrelevant for us because it will be cancelled in future 
computations. 

Then we define the partial determinants Di replacing the i—th column by the 
column of fi, putting them on the place of the n-th column: 



l<a<n 



-1 



%iz2) e,iz2) ... 

h{Zn) Oi{Zn) ... 

eo{z,) e,{z,) 
eo{z2) e,{z2) 



dn^lizi) /i 
On-l{z2) f2 

^n—li^Zji) fn 

■ On-lizi) 

■ On-l{z2) 



fa- 



Oo{Zn) Oi{Zri) ... On-l{Zn) 

Here the subscription | means that we had to omit the i-th column and the 
a- th row. 

The immediate corollary of the (1.1) is the following 

Proposition 3.1. The determinant ratios are formed a commutative family: 

[D-'D,,D-'Dj]^0. 

The result of the proposition can be expressed in an elegant way in the form 
of a commutation relation of generating functions. 

Let us define a generating function T{u) of a variable u e C: 

l<j<n 

Then we can express the function T(u), using the formulas for the determinants 
of the theta-functions as 



l<a<n 



eo{zi) ei{zi) 



On-l{zi) 



On-l{u) 



Oo{Zn) Oi{Zn) ... 9n-l{Zn) 



fa 



(12) 



e{:a+Y.^p) n O{u-Z0) 

E(3^a l<f^j^a<n 
fTluZ T^; 

l<a<n 



n o{za- Zfj) 

p^a 
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where we denote by fa the normahzation 

r fa 

ja n 

i=l 

We remark that the commutation relations between the variables Zi, . . . ,Zn,fi,--- , f, 
are the same as they were for the variables Zi, . . . ,Zn,fi,--- , fn- 
In this notations the proposition now reads: 

Proposition 3.2. The "transfer-like" operators T{u) commute for different val- 
ues of the parameter u: 

[TH,T(t;)]=0. 



Now we will apply this result to a construction of the commuting elements in 
the algebra Qn{^', v) even n. 

Let n = 2m. It is known in this case that the centre Z{Q2m{^',v)) for rj 
of infinite order is generated by Casimir elements from 5'™(62m(r)) such that 
f{zi, . . . , Zm) = for Z2 = Zi-\- 2mrj. A straightforward computation shows that 
the space of such elements is two-dimensional and has as a basis the elements of 
the form 



9a{zi + ■ ■ ■ ^m) JJ 0{Zi - Zj - 2m7]), 

where 6^ G 62(1), a G Z/2Z. 

Let us fix an element 'if{z) G Qm+^iX) and two complex numbers a, 6 G C. 
Consider a family f{u){zi, .... Zn) of elements from <S'™(02m(r)) such that 

f{u){zi, zi + 2mr], Z2 . . . , Zm-i) 

= "i/izi+im^T] ^{a+{m-2)b+2m{m-2)r]))e{zi+. . .+Zm-i+a)9{zi+Z2+b) . . . 

m + 5 

9{zi+Zm-i+b)9{z2—zi—4:mri) . . .9{zm-i—zi—Amr})9{z2—zi^-2mr]) . . .9{zm-i—zi+2mr})x 
9{u + Z2 + . . . + Zm-i — a — b + 2mri)9{u — Z2) ■ ■ ■ 9{u — Zm-i) x 
exp(27ri(2(m — 2)zi + Z2 + ■ ■ ■ + Zm-i)) Y\_ ^(^« ~ ~ 2mrf). 

2<i^j<m-l 

The elements f{u) are defined up to a linear combination of the Casimirs Ci, C2 
because of the annihilation of the Casimirs on the "diagonal" Zi — Z2 + 2mr] ( 
and we are defying the elements f{u) namely on this "diagonal"!) 

Theorem 3.1. In the elliptic algebra Q2mi^,v) the following relation holds 
[f(u), f(v)] = f(u) * f(v) - f(v) * f(u) = 
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Proof We will use the homomorphism Q2mi^,v) ~^ -Bm-i,2m from the sub- 
section 2.2.3. The element f{u) is transposed by this homomorphism into the 
element 

0{u+^Zf,) n 0{u-Zf,) 

,<4" 1 n o{zfs - 

l<a<m— 1 „, ^ ' 

- - P^a 

where we had denote by the the following expression 

= ^(2;^ + 4m^77 + (m - 2)6 + 2m(m - 2)77)) x 

m— 1 

^(Zl -^/^ + a) JJ^(^a + -2/3 + &) X 
13=1 fijta 

exp(27ri(2(m - 2)2:0 + ^ ^/3))eaei . . . e„_i. 

Then, the formula (12) and the proposition 3.2 give us immediately that the 

images of f{u) under the homomorphism (^^-i are commuting in Bm-i,2m{v) ■ It 
is known that the image of the algebra Q2m{S,v) -Bm-i,2m(^7) is the quotient 
of Q2m{^,v) over the centre Z = (Ci,C2). Hence the commutator 

[f{u)J{v)] = f{u)*f{v)-f{v)*f{u) 

belongs to the ideal generated by the centre Z. To show that [f{u),f{v)] = 
we can consider the injective homomorphism into Bm,2m{v) ^-nd it is sufficient 
to verify that the coefficient before (ei)^ . . . (e^)^ equals to zero, or ( which is 
equivalent) that 

[/(m), f{v)]{zi,zi + 2mr],Z2,Z2 + 2mr], ... ,Zm,Zm + 2mr]) = 0, 

which is achieved easily by the direct verification. 



Remark 3. A family of compatible quadratic Poisson structures such that its 
common linear combination is isomorphic to the Poisson structure of the classical 
elliptic algebra qn{S) was constructed in [31]. The Lenard scheme enables us 
with a family of Poisson commuting elements (" hamiltonians" in involution) in 
the algebra qn{£)- These elements have the degree n if n is impair and n/2 
otherwise. 

Let n = 2m. We conjecture that the constructed in 3.1 commuting elements 
in Q2m{£) arc the quantum analogs of the commuting "hamiltonians" in q2m{£) 
generated by the Lenard scheme applied to the classical Casimirs Co'-™-', Ci*^™^ of 
[31]. This conjecture is true in first non-trivial case n = 6 when it is easily verified 
by direct computations. 
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3.2. Commuting elements in "bosonization" of Qn,n-i{^,v)- The elliptic 
algebra Qn,n-i{^,v) is commutative and the homomorphism of "bosonization" 
provides a construction of a large class of the commuting families in the corre- 
sponding Weyl-like algebras. Some special choice of the relations between the 
numeric parameters (numbers of the generators) gives us in its turn the example 
of some new (to our knowledge) integrable system. 

3.2.1. Bosonization of the algebra Qn n_i{S) . . The homomorphism 0^ : Qn{S,ri) - 
Bn,p{i]) from 2.2.3 may be extended and generalized to the case of the elliptic al- 
gebras Qn,k{^,ri) ([27]). The structure of the Weyl-like algebra similar to Bn^p{rj) 
turns out to be more complicated for /c > 1. We will use this extension in the 
special case k = n — 1 i.e. for the commutative algebra Qn,n-i{£, 

Let us consider the expansion of in continued fraction to obtain the numeric 
parameters of the dynamical Weyl-like algebra similar to Bn,p{r]): 



n-1 2-...-\- 

Let Bn^p^^p^^,,. ^p^_-^{r]) be an associative algebra generated by the following gen- 
erators 

-2^1,1) • • • ) • • • ) 6pi,lj 

^1,2, ■ ■ • , 2:^2,2; ei,2, • • • , ep2,2; : 

^l.n— 1) • • • ; Zp„-i,n—li 6l,n— 1) ■ ■ • ) ^pn-i,n—li 

and the additional generators 

^1,2, ^2,3 • • • , tn-2,n-l', 
/l,2, /2,3 • • • • fu-2,n-l', 

We will impose the following commutation relations between them: 

{ta,a+l - 'nr))fa,a+l- 

All other commutators we suppose to be equal to zero. 

The commutative elliptic algebra Qn,n-ii^,i]) admits the following "bosoniza- 
tion" homomorphism 0pi,...,p„_i which is a partial case of the general bosonization 
procedure for Qn,k{^,v)- 

More precisely, 

can be squeezed through a homomorphism of the algebra Qn-i,n into an alge- 
bra of "exchange relations" with generators eci,.. (see [30]). This "exchange" 
algebra, in its turn, can be imbed in the algebra -Bpi,...,p„_i(?7) with additional 
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generators. This statement was established in the paper [28]. The exphcit com- 
position map (ppi,... ,pn-i t)e expressed Uke the following "transfer-function" 

)e{u + )...0{u + Za„_^,n-l) 

l<a,<p, n ^(^a„7 - ^At) (13) 

- - /3 # a-, ^ ^ 

1 < 7 < n — 1 

+ ^Q2,2 — ^1,2)^(2:02,2 + 2:0,3,3 — t2,3) • • • ^(2:a„_2,„-2 + Zan-i,«-i " ^n-2,n-l) 
6ai,l ■ ■ • Ca„_i „_i/l,2 ■ ■ ■ fn-2,n-l- 

Then the following proposition is resulted from the definition of the combined 
homomorphism of the bosonization: 

Proposition 3.3. 

[f{u),f{v)]^0. 

We should observe that in the case of pi = . . . = Pn-i — 2 these commuting 
family gives an example of integrable system. 

4. Some examples of elliptic integrable systems 

4.1. Low- dimensional example: Algebra q2{£,r])- The elliptic integrable 
system arising in the 3.1 becomes very transparent in the case of m = 1 then the 
commutative elliptic algebra Q2{£,t]) has functional dimension 2 and its Poisson 
counterpart admits the Poisson morphism 

i)2 ■ q2{S) 62, 

where the algebra 62 are formed by the elements 

X]/a,/3(^i,^2)eief, 

a,/3 

where fa,i3{zi, Z2) are meromorphic functions and the Poisson structure in 62 is 
given by 

{zi,Zj} = {ei,ej} = {ei,Zi} = 0{ei,Zj} = -2ei{i = 1,2. (14) 

The explicit formula for this mapping (for a given theta-function of the order 2 
/ e 62 (F)) is the following: 

/ r f\ /(^l) I /(^2) 

M) = 77; -rei + -^62- (15) 

&[Zi - Z2) U[Z2 - Zi) 

Now, taking the basic theta-functions of the order 2 61,62, let us compute the 
Poisson brackets between their images ■02 (^1) and 1^2(62) ■ 

Proposition 4.1. These theta-functions commute in 62 
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4.2. SOS eight-vertex model of Date-Miwa-Jimbo-Okado. Wc remind in 
a neccessairy form some of tlie SOS eight-vertex model's ingredients (see [23] ) 
and the construction of its transfer-operators to establish the relation between 
them and our operators T{u). 

This is a statistical mechanical model which is an Interaction-round- a- face - 
model (or IRF-model )-version of the Baxter model related to the elliptic quantum 
group of Felder ET^nish) which was studied by the Sklyanin's Separartion of 
Variables methods (under the antiperiodic boundary conditions) in [25] (see also 
[24] for the representation theory of £'^^^(5/2))- We will use the results of [25] in 
the form which we need. 

The antiperiodic boundary conditions of the model are fixed by the family of 
transfer-matrix Tsosiu, A) where w e C is a parameter and the family Tsos{u, A) 
is expressed as a (twisted) traces of the L— operators of the model Lsosiu, A) de- 
fined over so-called "auxiliary" space of the elliptic quantum group Er^-qish)- The 
L— operator is built out in the tensor product of the fundamental representations 

of the elliptic quantum group and is twisted by the matrix 

The L— operator is usually represented in the form of 2 x 2 — matrix of the 
form 




Lsos{u, A) 



a{u, A) b{u, A) 
c{u, A) d{u, A) 



where the matrix entries are meromorphic on u and A and obey the dynamical 
it!LL— commutation relations for the Felder elliptic i?— matrix 



R{u,X) 



( e{u + 277) 



V 









e{u)e(\+'2ri) 


d{u-\)d{2r]) 


e{\) 


e(x) 




0{u)e(X-2ri) 



















9{u + 29) J 



(see op.cit). 

We will have deal with the functional representations of £'^-^^(5/2) which are 
the pairs {F, L) where F is a complex vector space of meromorphic functions 
f{zi, ... ,Zn,X) ( or a subspace of functions which are holomorphic in a part of the 
variables) and L is the L— operator as above. The entries of the L— operator are 
acting as the difference operators in the tensor product V <SiW where dimV — 2 
such that the Felder i?— matrix belongs to End(y (8) V) and W is an appropriate 
subspace in the functional space F. 

The Bethe ansatz method works for case of SOS model with periodic boundary 
conditions after Felder and Varchenko ([24]). In the antiperiodic case the family 
of the transfer matrices 



Tsos{u, A) = trKLsos{u, A), ueC 
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is commutative for the different parameter values: [Tsosi^, X) jTsosi^, X)] = 0, 
as it follows from the RLL— relations by tediuos computations in [26]. 

In the other hand it is possible to estabhsh the exphcit one-to-one correspon- 
dence between the families of antiperiodic SOS transfer-matrices and the auxiliary 
transfer-matrices Taux{u, \) (see 4.4.3 in [26]). This isomorphism is established 
by the version of the Separation of Variables. 

The explicite expression of the auxiliar transfer-matrix is 

where ( to compare with our "transfer" -operators in the elliptic intgerable sys- 
tems) we had put in the formulas of ([26], ch.4.6) 

Xa^O, Aa = 1, a; = 1, ...,n. 

(The choice = 1 is corresponded to the case Separated Variables (prop. 4.36 
in [26])). 

The operators are acting in the way similar to the generators above: 

^n) = f{zi, ± 277, Zn)T^^\ 
Let us make the following change of the variables 

Now the simple inspection of the formulas shows that 

Tauxiu,X)^T+{u,X)+T4u,X), 
and the the commutation results of 3.2 can be applied by the following 

Proposition 4.2. The " transfer"- operator T{u) coincides (up to an inessential 
numeric factor, depending on the normalization in the theta-function definition 
and rescaling of the parameter r]) with the comMnantion of the twisted traces 
of the auxiliary L— operator of the antiperiodic SOS-model under the following 
choice of the variables: 

n 
a=l 

The proposition 3.2 gives a simple proof of the commutation relations for the 
twisted traces of the SOS-model in this special case. We should observe that the 
commutation relations 

[T4u),T4v)] + [T4u),T4v)]^0 
are followed from the same arguments like in 3.2. 
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In other hand, we have got an additional argument to our behef that the 
"integrabihty" condition from [12] is equivalent in some sense to the another 
"integrability" form ( in the given case to the i?LL— relations). 

We beheve also that the role of the elliptic integrable systems in the IRF models 
does not restrict only to some specific conditions and we hope to return to this 
point in future. 

4.3. Elliptic analogs of the Beauville-Mukai systems and Fay identity. 

Let us describe a geometric meaning of our IS's. We will start with an observation 
that the classical analog of the Weyl-like algebra Vn given by the Poisson brackets 

= {fh fj} = {zi, Zj} = {fi, Zj} {i j), {fi, Zi} = -nfi 
can be identified with the Poisson algebra of the functions on the symmetric prod- 
uct S"'{Cone{S)) of a surface represented the elliptic curve cone ( more precisely, 
on the Hilbert scheme {Cone{S))^'^^ of the length n of the points on this surface.) 

The n commuting elements hi — D^^Di in the Poisson algebra obtaining in 
the 3.2 can be interpreted as an elliptic version of the Beauville-Mukai systems 
associated with this Poisson surface. 

We will develop in more details the first interesting case (n = 3) of the Poisson 
commuting conditions (the Plucker relations from sect. 2) for this systems. 

The Beauville-Mukai hamiltonians have the following form in this case: 

det[e,^i,^2] rr det[e,^o,^2] „ det[e,^o,^i] 



det[^o,^i,^2]' det[^o,^i,^2]' det[^o, ^i, ^2] ' 
where the vectors-columns e,9i, i — 0,1,2 have the following entries 



e^= ft.(^2) , i = 0,l,2 





Now the integrability condition (6) may be exprimed as a kind of 4-Riemann 
identity (known also as the trisecant Fay's identity): 

emu . /> . j% ^ . jUsiu . ^^^^^ 

'■^ E{A,D)E{C,B) 



0{u+ / uj)e{u+ / cu) 

JB Jd 



E{A, C)E{D,By 



where : 8 ® £ ^ Jac(£) = £ is a Jacobi mapping X,Y ^ 8, := uj 
for a chosen holomorphic differential uj on the elliptic curve, 
E{X, Y) = -E{Y, X) is a prime form: 

E{A,B)E{D,c) _~e{!l)0U^) 



E{A,c)E{D,B) e{j^)e{j;^) 
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and 6{u) is an odd thcta- function, (sec [45].) 

Let a, b, c, d be four points corresponded to the point A, B, C, D in (16) and we 
will choose u — a + c. 
Then the (16) reads as 

e(a + c)e{a - c)e{h + d)e{h -d)- e{a + b)e{a - b)e{c + d)e{c -d) + 

(17) 

e{a + d)e{a - d)e{c + b)e{c - &) = 0, 

where we are easily recognizing a form of the commutativity conditions (6) for 

the case n = 3 and an appropriate choice of a linear relation between 4 points 
a, 6, c, d, and Zi^ Z2. z^.t] /?> (modulo an irrelevant exponential factor entering in the 
relations between the theta-functions in different normalizations). 

Remark 4. 1. The appearance of the Riemann-Fay relations as a commutativ- 
ity or an integrability conditions in this context looks quite natural both from 
the "Poissonian " as well as from the "integrable" viewpoints. The "Poisson- 
Plucker" relations and their generalizations in the context of the Poisson polyno- 
mial structures were studied in [44]. In the other hand, the first links between the 
integrability conditions ( in the form of Hirota bilinear identities for some elliptic 
difference many-body-like systems) and the Fay formulas were established in [17]. 

2. Second remark concerns the relations between the Fay trisecant formulas 
on an elliptic curve and a version of a "triangle" relations known as the "asso- 
ciative" Yang-Baxter equation obtained by Polischuk ([42]). This result gives an 
additional evidence that the commutation relations (2), (3) could be interpreted 
as an algebraic sort of Yang-Baxter equation. An amusing appearance of the 
NC determinants in both constructions (the Cartier-Foata determinants defined 
above are of course a partial case of the quasi-determinants of Gelfand-Retakh [7]) 
shows that the ideas of [12] might be useful in "non-commutative integrability" 
constructions which involve quasi-determinants, quasi-Plucker relations etc. 

5. Discussion and future problems 

We have proved that the Elliptic Algebras ( under some mild restrictions) carry 
the families of commuting elements which become in some cases the examples of 
Integrable Systems. 

Let us indicate some questions deserving future investigations. 

We are going to construct an analogs of the IS on the algebras Q2m{^,v) 
which are corresponded to the maximal symplectic leaves of the Elliptic Alge- 
bras Q2m+i{^,v) using the bi-hamiltonian elliptic families. This analog should 
quantize the bi-hamiltonian families in C" in the case of an impair n. 

One of the main motivations to construct the systems on Vn explicitly in the 
terms of the determinants of functions of order n was the tentative to find 
a confirmation to the hypothetical "separated" form of the hamiltoninans to 
n— point Double- Elliptic system, proposed in the paper [4] . 
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The relevance of our construction to this circle of the problems had got an ad- 
ditional evidence in the paper [35] . Recent discussions around of the integrability 
in Dijkgraaf-Vafa and Seiberg-Witten theory provide us with the "physical" in- 
sights to support the relation between the Beauville-Mukai and Double Elliptic 

IS ( see [6], [5]). 

The future paper (in collaboration with A. Gorsky) which should clarify the 
place of the IS associated with the Elliptic Algebras and the integrability phe- 
nomena in SUSY gauge theories are in progress. 

There are some open questions about the relations of these IS with the Beauville- 
Mukai Lagrangian fibrations on the Hilbert scheme (P^yi")'"! as well as with their 
non-commutative counterparts proposed in [43]. We can argue that the proposed 
in [9] quantization scheme for the fraction fields can be apphed to the NC surface 
(P5\£) and to "extend" the deformation to the NC Hilbert scheme (P^y^")!"' in- 
troducing in [43]. The proposition 3.1 from [12] then should in principle to give 
an NC integrable Beauville-Mukai system on (P^y^^)'"' as a NC IS associated 
with the Cherednik algebras of [10]. These questions are also in the scope of our 
interest. 

Finally, we should mention an interesting question of generalization of the 
Beauville-Mukai IS which was introduced by Gelfand-Zakharevich on the Hilbert 
scheme (Xg) f"l of the Del Pezzo surface Xg getting by the blow-up of 9 points on 
P^. These 9 points are the intersection points of two cubic plane curves. Let tt is 
a Poisson tensor on Xg then it can be extend to the whole P^ such the extension 
TT has 9 zeroes in these 9 points. The polynomial degree of the tensor tt is 3 and 
it follows that this polynomial is a linear combination of the polynomials corre- 
sponding to the initial eUiptic curves. This bi-hamiltonian structure corresponds 
to the case n — 3 studied in the paper [31]. The open question arises - what 
an algebro-geometric construction behind the bi-hamiltonian Elliptic Poisson Al- 
gbera of [31] in the case of arbitrary nl 
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